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We present measurements of branching fractions and CP-violating asymmetries in B° — ► p ± n^ i 
and B° — > p~ K + decays. The results are obtained from a data sample of 88.9 x 10 6 T(AS) — > BB 
decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. 
From a time-dependent maximum likelihood fit we measure the charge-averaged branching fractions 
B(B° -> p ± Tx^) = (22.6±1.8 (stat)±2.2(syst))x!0- 6 and B(B a -> p~K+) = (7.3 ti'j ± 1.3) x 10" 6 ; 
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and the CP-violating charge asymmetries A p J p = -0.18 ± 0.08 ± 0.03 and A p c K p = 0.28 ± 0.17 ± 0.08; 
the direct CP violation parameter C p7T = 0.36 ± 0.18 ± 0.04 and the mixing- induced CP violation 
parameter S pn = 0.19 ± 0.24 ± 0.03; and the dilution parameters AC pn = 0.28± o ;ig ± 0.04 and 
AS P n =0.15 ±0.25 ±0.03. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



In the Standard Model, CP-violating effects arise from 
a single complex phase in the three-generation Cabibbo- 
Kobayashi-Maskawa quark- mixing matrix 1] . One of the 
central, unresolved questions is whether this mechanism 
is sufficient to explain the pattern of CP violation ob- 
served in nature. We present here a simultaneous mea- 
surement of branching fractions and CP- violating asym- 
metries in the decays B° — > /9 ± 7r =F and B° — ► p~K + 
(and their charge conjugates). The BABAR and Belle 
experiments have performed searches for CP-violating 
asymmetries in B decays to 7r + 7r~ 0, 0, where the 
mixing-induced CP asymmetry is related to the angle 
a = Big[-V td V^/V ud V* b ] of the Unitarity Triangle as 
it is for p^tt^. However, unlike 7r + 7r~, p^n^ is not 
a CP eigenstate, and four flavor-charge configurations 
(B°(B°) — > p ± 7r =F ) must be considered. Although this 
leads to a more complicated analysis, it benefits from a 
branching fraction that is nearly five times larger 0, 0] . 

Following a quasi- two-body approach |2j, we restrict 
the analysis to the two regions of the n^^h^ Dalitz 
plot (h = 7T or K) that are dominated by either p + h~ 
or p~ h + . With At = t p h — t tag defined as the proper time 
interval between the decay of the reconstructed B° h and 
that of the other meson P°agj the time-dependent decay 
rates are given by 



■\At\/r 
~4t 



(1) 



1 ± Qt m (S ph ± AS ph ) sm{Am d At) 
- Qtag (C ph ± AC ph ) cos(Am d At) 



3° 

5 tae 



IS 



where Qtag = 1( — 1) when the tagging meson B. 
a B°(B°), t is the mean B° lifetime, and Arrid is 
the B°B° oscillation frequency. The time- and flavor- 
integrated charge asymmetries A^ p and A p ^ p measure 
direct CP violation. For the pn mode, the quanti- 
ties S plT and C p7T parameterize mixing-induced CP vi- 
olation related to the angle a, and flavor-dependent di- 
rect CP violation, respectively. The parameters AC p7r 
and AS p7T are insensitive to CP violation. AC p7T de- 
scribes the asymmetry between the rates r(P° — > 
p+TT-) + T{B° -> p-ir+) and r(P° -> p-ir+) + T(B° -> 
p + 7T _ ), while AS p7r is related to the strong phase dif- 
ference between the amplitudes contributing to B° — > 
pn decays. More precisely, one finds the relations 
S /m ±AS /m = y/1- (C p7r ±AC p7r ) 2 sin(2c4 ±5), where 
2a± = arg[(<?/p)(A± /A*)], S = axg[A-/A+], argfe/p] 



is the B°B° mixing phase, and A+ (A+ ) and A p7r (A p7I ) 
are the transition amplitudes of the processes B°(B°) — > 
p + 7r~ and B°(B°) — > p~n + ', respectively. The angles cv^ 
are equal to a in the absence of contributions from pen- 
guin amplitudes. For the self-tagging mode, the val- 
ues of the four time-dependent parameters are C p k = 0, 
AC p k = -1, Sp K = 0, and AS pA - = 0. 

The data used in this analysis were accumulated with 
the BABAR detector ||, at the PEP-II asymmetric- 
energy e + e~ storage ring at SLAC. The sample consists 
of (88.9 ±1.0) x 10 6 BB pairs collected at the Y(4S) 
resonance ( "on- resonance" ) , and an integrated luminos- 
ity of 9.6 ftT 1 collected about 40 MeV below the Y(4S) 
( "off- resonance" ) . In Ref. Q we describe the silicon ver- 
tex tracker and drift chamber used for track and vertex 
reconstruction, the Cherenkov detector (DIRC) , the elec- 
tromagnetic calorimeter (EMC), and their use in particle 
identification (PID). 

We reconstruct B ph candidates from combinations of 
two tracks and a ir° candidate. We require that the PID 
of both tracks be inconsistent with the electron hypoth- 
esis, and the PID of the track used to form the p be 
inconsistent with the kaon hypothesis. The ir° candidate 
mass must satisfy 0.11 < 771(77) < 0.16 GeV/c 2 , where 
each photon is required to have an energy greater than 
50 MeV in the laboratory frame and to exhibit a lateral 
profile of energy deposition in the EMC consistent with 
an electromagnetic shower. The mass of the p candidate 
must satisfy 0.4 < r7i(7r ± 7r°) < 1.3 GeV/c 2 . To avoid the 
interference region, the B candidate is rejected if both 
the 7r + 7T° and 7r _ 7r° pairs satisfy this requirement. Tak- 
ing advantage of the helicity structure of B — > ph de- 
cays (h is denoted bachelor track hereafter), we require 
| cos 6*^1 > 0.25, where 6„ is the angle between the 7r° 
momentum and the negative B momentum in the p rest 
frame. The bachelor track from the ph decay must have a 
e + e~ center-of-mass (CM) momentum above 2.4GeV/c. 

For 86% of the B a — > ph decays that pass the event 
selection, the pion from the p has momentum below this 
value, and thus the charge of the p is determined unam- 
biguously. For the remaining events, the charge of the p 
is taken to be that of the tt^tt combination with mass 
closer to the p mass. With this procedure, 5% of the se- 
lected simulated signal events are assigned an incorrect 
charge. 

To reject background from two-body B decays, the in- 
variant masses of the n^h^ and h^ir combinations must 
each be less than 5.14 GeV/c 2 . Two kinematic variables 
allow the discrimination of signal- P» decays from fake-B 



5 



decays due to random combinations of tracks and it can- 
didates. One variable is the difference, AE, between the 
CM energy of the B candidate and y/s/2, where y/s is the 
total CM energy. The other variable is the beam-energy- 
substituted mass m,Es = \J (s/2 + Pi ■ ps) 2 /Ef — p 2 B , 
where the B momentum ps and the four- momentum of 
the initial state (Ei, p,) are defined in the laboratory 
frame. The AE distribution for pir (pK) signal peaks 
around (—45) MeV since the pion mass is always as- 
signed to the bachelor track. We require 5.23 < tues < 
5.29 GeV/c 2 and -0.12 < AE < 0.15 GeV, where the 
asymmetric AE window suppresses higher-multiplicity B 
background, which leads to mostly negative AE values. 
Discrimination between pn and pK events is provided by 
the Cherenkov angle &c and, to a lesser extent, by AE. 

Continuum e + e~ — > qq (q = u,d,s,c) events are the 
dominant background. To enhance discrimination be- 
tween signal and continuum, we use a neural network 
(NN) to combine four discriminating variables: the re- 
constructed p mass, cos 6^, and the two event-shape vari- 
ables that are used in the Fisher discriminant of Ref. Q . 
The NN is trained in the signal region with off-resonance 
data and simulated signal events. The final sample of sig- 
nal candidates is selected with a cut on the NN output 
that retains ~ 65% (5%) of the signal (continuum). 

Approximately 23% (20%) of simulated pir (pK) events 
have more than one ph candidate passing the selection 
criteria. In these cases, we choose the candidate with the 
reconstructed ir° mass closest to the nominal ir° mass. 
A total of 20,497 events pass all selection criteria. The 
signal efficiency determined from Monte Carlo (MC) sim- 
ulation is 20.7% (18.5%) for pir {pK) events; 31% (30%) 
of the selected events are misreconstructed, mostly due 
to combinatorial- 7T° background. 

We use MC-simulated events to study the cross-feed 
from other B decays. The charmless modes are grouped 
into eleven classes with similar kinematic and topolog- 
ical properties. Two additional classes account for the 
neutral and charged b — > c decays. For each of the back- 
ground classes, a component is introduced into the like- 
lihood, with a fixed number of events. In the selected pir 
(pK) samples we expect 6±1 (20±2) charmless two-body 
background events, 93±23 (87±22) charmless three-body 
background events, 118±65 (36±18) charmless four-body 
background events, and 266 ± 43 (54 ± 11) b — ► c events. 
Backgrounds from two-, three-, and four-body decays to 
pn are dominated by B + — > 7r + 7r°, B + — > p°tt + , and 
longitudinally polarized B° — > p + p~ decays. The pK 
sample receives its dominant two-body background from 
B + — > K + ir°, and its dominant three- and four-body 
background from B — > K*tt and higher kaon resonances, 
estimated from inclusive B — + Ktttt measurements. 

The time difference At is obtained from the measured 
distance between the z positions (along the beam direc- 
tion) of the B° h and B® ag decay vertices, and the boost 
/?7 = 0.56 of the e + e~ system 0,0- To determine the 



flavor of the B® ag we use the tagging algorithm of Ref. . 
This produces four mutually exclusive tagging categories. 
We also retain untagged events in a fifth category to im- 
prove the efficiency of the signal selection and the sensi- 
tivity to charge asymmetries. Correlations between the 
B flavor tag and the charge of the reconstructed ph candi- 
date are observed in various i?-background channels and 
evaluated with MC simulation. 

We use an unbinned extended maximum likelihood fit 
to extract the pn and pK event yields, the CP param- 
eters and the other parameters defined in Eq. (JJJ. The 
likelihood for the N% candidates i tagged in category k 
is 

N k tt,K , N B 

c * = e Nk Ii E NPh ^t k + N?*P% h + £ 

i=l h j=\ ' 

where N' k is the sum of the signal and continuum yields 
(to be determined by the fit) and the fixed B-background 
yields, N ph is the number of signal events of type ph 
in the entire sample, tk is the fraction of signal events 
tagged in category k, and N^' h is the number of con- 
tinuum background events with bachelor track of type h 
that are tagged in category k. The total likelihood C is 
the product of likelihoods for each tagging category. 

The probability density functions (PDFs) V p k h , T>f ,h 
and the likelihood terms £f are the product of the 
PDFs of five discriminating variables. The signal PDF 
is thus given by V p k = V ph (m ES ) • V ph (AE) ■ V ph (NN) ■ 
V ph (6c)-'Pu (At), where 'P k h ( At) contains the measured 
physics quantities defined in Eq. l(T|l. diluted by the ef- 
fects of mistagging and the At resolution. The PDF of 
the continuum contribution with bachelor track h is de- 
noted 'P k ™ ,h . The likelihood term L® k h corresponds to 
the i?-background contribution j of the Nb categories. 

The signal PDFs are decomposed into three parts with 
distinct distributions: signal events that are correctly re- 
constructed, misreconstructed signal events with right- 
sign p charge, and misreconstructed signal events with 
wrong-sign p charge. Their individual fractions are esti- 
mated by MC simulation. The rriES, AE, and NN output 
PDFs for signal and B background are taken from the 
simulation except for the means of the signal Gaussian 
PDFs for ities and AE, which are free to vary in the fit. 
The continuum PDFs are described by six free param- 
eters. The B c PDF is modeled as in Ref. 0|. The Af- 
resolution function for signal and B-background events 
is a sum of three Gaussian distributions, with parame- 
ters determined from a fit to fully reconstructed B° de- 
cays 0] . The continuum At distribution is parameterized 
as the sum of three Gaussian distributions with common 
mean, two relative fractions, and three distinct widths 
that scale the At event-by-event error, yielding six free 
parameters. For continuum, two charge asymmetries and 
the ten parameters N k ^ ,h are free. A total of 34 param- 
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TABLE I: Summary of the systematic uncertainties. 



Error source 


(events) 


AP K 

A CP 


A CP 

(in 


C p-K AC p-K 
units of 10 


) 




Amd and r 


0.1 


0.1 


0.0 


0.0 


0.4 


0.4 


0.2 


0.1 


At PDF 


1.2 


1.9 


0.4 


0.2 


1.4 


0.8 


1.5 


1.2 


Signal model 


4.0 


13.1 


1.2 


0.8 


0.7 


0.8 


1.4 


1.0 


Particle ID 


0.6 


0.7 


0.5 


0.2 


0.1 


0.1 


0.1 


0.1 


Fit procedure 


8.0 


15.7 


0.4 


0.2 


0.4 


0.4 


0.4 


0.3 


DCS decays 


0.0 


0.3 


0.0 


0.1 


2.2 


2.2 


0.8 


0.7 


B background 


16.0 


14.2 


7.9 


2.8 


3.0 


3.5 


2.1 


1.8 


Total 


18.4 


25.0 


8.0 


2.9 


4.1 


4.3 


3.1 


2.5 



eters, including signal yields and the parameters from 
Eq. 0, are varied in the fit. 

The contributions to the systematic error on the signal 
parameters are summarized in Table [i] The uncertain- 
ties associated with Arrid and r are estimated by vary- 
ing these parameters within the uncertainty on the world 
average [lfj. The uncertainties due to the signal model 
are obtained from a control sample of fully reconstructed 
B° — ► D~ p + decays. We perform fits on large MC sam- 
ples with the measured proportions of pit / pK signal, and 
continuum and B backgrounds. Biases observed in these 
tests are due to imperfections in the PDF model; e.g., un- 
accounted correlations between the discriminating vari- 
ables of the signal and B-background PDFs. The biases 
are added in quadrature and assigned as a systematic un- 
certainty of the fit procedure. The systematic errors due 
to interference between the doubly-Cabibbo-suppressed 
(DCS) b — * ucd amplitude with the Cabibbo-favored 
b — > cud amplitude for tag-side B decays have been es- 
timated from simulation by varying freely all relevant 
strong phases [llj . 

The main source of systematic uncertainty is the B- 
background model. The expected event yields from the 
background modes are varied according to the uncer- 
tainties in the measured or estimated branching frac- 
tions. Systematic errors due to possible nonresonant 
B° —> ir + ir~ir° decays are derived from experimental lim- 
its Repeating the fit without using the p-candidate 
mass and helicity angle gives results that are compatible 
with those reported here. Since £?-background modes 
may exhibit CP violation, the corresponding parameters 
are varied within their physical ranges. 

The maximum likelihood fit results in the event yields 
N p * = 428 tH and N pK = 120 til, where the errors 
are statistical. Correcting the yields by a small fit bias 
determined using the MC simulation (3% for pn and 0% 
for pK), we find for the branching fractions 

B(B° ^ p ± n T ) = (22.6 ± 1.8 ± 2.2) x 1CT 6 , 
B{B"^p-K + ) = (7.3±H± 1.3) x 1CT 6 , 

where the first errors are statistical and the second sys- 




m „ GeV/c rn „ GeV/c 




-0.1 -0.05 0.05 0.1 0.15 -0.1 -0.05 0.05 0.1 0.15 



AE GeV AE GeV 

FIG. 1: Distributions of tues and AE for samples enhanced in 
pn signal (a,c) and pK signal (b,d). The solid curve represents 
a projection of the maximum likelihood fit result. The dashed 
curve represents the contribution from continuum events, and 
the dotted line indicates the combined contributions from con- 
tinuum events and B-related backgrounds. 

tematic. The systematic errors include an uncertainty of 
7.7% for efficiency corrections, dominated by the uncer- 
tainty in the 7r° reconstruction efficiency. Figure ^ shows 
distributions of tues and AE, enhanced in signal content 
by cuts on the signal-to-continuum likelihood ratios of 
the other discriminating variables. For the CP- violating 
parameters, we obtain 

A P J P = -0.18 ± 0.08 ± 0.03 , A P C K P = 0.28 ± 0.17 ± 0.08 , 
C l „= 0.36 ±0.18 ±0.04, Spn =0.19 ±0.24 ±0.03. 

For the other parameters in the description of the 
B°(B°) — > pir decay-time dependence, we find 

ACV = 0.28±2;i| ± 0.04 , AS p7r = 0.15 ± 0.25 ± 0.03 . 

We find the linear correlation coefficients cc,ac = 0.18 
and cs,as — 0.23, while all other correlations are smaller. 
As a validation of our treatment of the time depen- 
dence we allow r and Ato<j to vary in the fit. We find 
t = (1.64 ± 0.13) ps and Am d = (0.52 ± 0.12) ps" 1 ; the 
remaining free parameters are consistent with the nom- 
inal fit. The raw time-dependent asymmetry A B o jpa — 
(N B o—Np (l )/(N B a +Ngo) in the tagging categories dom- 
inated by kaons and leptons is represented in Fig. [21 

In summary, we have presented measurements of 
branching fractions and CP-violating asymmetries in 
B° — > p ± iT T and p~K + decays. We do not observe direct 
or mixing-induced CP violation in the time-dependent 
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FIG. 2: Time distributions for events selected to enhance 
the p7r signal tagged as (a) -B? ag and (b) -E>t ag , an< i ( c ) time- 
dependent asymmetry between -Bta g and -B?a g - The solid 
curve is a likelihood projection of the fit result. The dashed 
line is the total B- and continuum-background contribution. 



asymmetry of B° — > p^n^ decays and there is no evi- 
dence for direct CP violation in B° — ► p~K + . 
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